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Recently several works have appeared in the literature in which authors try to describe 
Freeze Out (FO) in energetic heavy ion collisions based on the Boltzmann Transport 

in 
o 



Equation (BTE). The aim of this work is to point out the limitations of the BTE, when 
applied for the modeling of FO or other very fast process, and to propose the way how 
p ■ the BTE approach can be generalized for such a processes. 

PS. 

The Freeze Out (FO) is an important phase of dynamical reactions. The FO is a kinetic 
process and thus the phase-space (PS) distribution of post FO particles can be obtained 
from kinetic FO calculations. The connection of the kinetic description of this process 
^ ■ and the Boltzmann Transport Equation (BTE) raised considerable attention recently 
|l|2j . One would think that FO can be handled perfectly by using the BTE, which may 
describe both equilibrium and non-equilibrium processes in a 4-dimensional space-time 
(ST) volume element like FO layer. Our aim is to analyse the situation and point out 
the physical causes, which limit the applicability of the BTE for describing FO and some 
other processes j^j. 

FO is usually assumed to happen on sharp 3-dimensional ST hypersurfaces. However, 
the FO-fronts or FO-layers are not necessarily narrow, but typically they have a charac- 
teristic direction or normal, <i 3 cr M , which can be time-like or space-like. It is more realistic 
to assume a continuous, 4-volume FO in a layer (or domain) of the ST. At the inside 
boundary of this layer there are only interacting particles, while at the outside boundary 
hypersurface all particles are frozen out and no interacting particles remain. Assuming 
that the boundaries of this layer are approximately parallel and neglecting derivatives in 
all directions other than g? 3 <t m , we can consider FO process as 1-dimensional |3|5j . 

We can derive the BTE from the conservation of charges in a ST domain, A 4 x, assuming 
the standard conditions, see for example [I]: 

(i) only binary collisions are considered, 

(ii) we assume "molecular chaos", i.e. that the number of binary collisions at position x 
is proportional to f(x,px) x /(x,p 2 ), 

(iii) f(x,p) is a smoothly varying function compared to the mean free path (m.f.p.). 



2 



We have to take into account that particles can scatter into the phase space (PS) volume 
element around p, or can scatter out from this volume element, described by Gain- and 
Loss- collision terms in the BTE. In these terms we consider elementary binary collisions 
where in the initial state two particles collide with momenta p\ and P2 into a final state of 
two particles with momenta p% and p^. In case of the Gain term the particle described by 
the BTE, with momentum p (without an index), is one from the two final state particles, 
while in case of the loss term this particle is one of the initial state particles. This is 
indicated by the indexes of the invariant transition rate, W 3 l and W 3 2 correspondingly 
[3]. We integrate over the momenta of the other three particles participating in this binary 
collision. We use the following notation for the PS integrals: 

_ d 3 p 1 d 3 p 2 d 3 p 3 

Urn = — o o o~ • 

Pi pi Ps 

So finally we can write BTE in the following way: 

P»d,f(x,p) = i J D 2U w£ [f{x,p 3 )f(x,p 4 ) - f(x,p)f(x,p 2 )} . (1) 

Here we want to underline the important symmetry properties of the invariant transition 
rate, namely Wf^ = = W^. This transition rate can be expressed as Wf^ = 

S(Tdiff8(Pi + P2 — P3 — Pi), where ^-function enforces energy-momentum conservation, 
s = (Pi + P2) 2 and Udiff is a differential cross section, here we will assume that it is 
momentum independent (so the total cross section a = Anadiff)- 

However, the usual structure of the collision terms in the BTE is not adequate for 
describing rapid process in a layer with a thickness comparable with the m.f.p. If we 
assume the existence of such a layer this immediately contradicts assumption (iii): the 
change is not negligible in the direction of d 3 a fl (normal to the layer). The assumption of 
"molecular chaos" (ii) is also violated in such a process, because number of collisions is 
not proportional to f(x,pi) x f(x,p 2 ), but it is delocalized in the normal direction with 
f(xi,pi) x f(x 2 ,p 2 ), where x^ is the origin of colliding particles, i.e., the ST point where 
the colliding particles were colliding last. 

Based on the above considerations, one might conclude that the changes of the distri- 
bution function are mediated by the transfer of particles, and consequently only slowly 
propagating changes are possible. If the FO layer propagates slowly, then its normal, 
d 3 a fl , must always be space-like. This was a common misconception, where all "superlu- 
minous" shock, detonation, deflagration fronts or discontinuities, and FO were considered 
unphysical based on early studies j^j. However, later it was shown that discontinuous 
changes may happen simultaneously in spatially neighbouring points, i.e. the normal of 
the discontinuity-hypersurface can be time-like [7]. This applies to the FO process also. 
Thus, the direction of characteristic or dominant change, d 3 a' M , may be both space-like 
and time-like in the FO process. 

From the all processes mentioned above (i.e. shocks, detonations, deflagrations etc.) 
the FO is the most special one, because the number of interacting particles is constantly 
decreasing as the FO proceeds, correspondingly the m.f.p. is increasing. In fact, it reaches 
infinity when the FO is completed. This means that we can not make the FO in finite 
layer of any thickness smooth enough to be modeled with the BTE! It is also obvious that 
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if FO has some characteristic length scale, it is not proportional to the m.f.p., because 
the m.f.p. increases as the density of interacting component becomes smaller, while the 
FO becomes faster in this limit, so its characteristic scale should decrease. 

To describe that the PS distributions which change rapidly, faster than the m.f.p., we 
suggest a Modified Boltzmann Transport Equation (MBTE): 

P»d,f(x, P ) = l -J d 2U wil x W~p^) x - W~p) x W~p~2) x ] , (2) 

where f(x,pi) x is an average over all possible origins of the particle in the backward 
lightcone of the ST point x = (t,x): 

where 5 3 (x — X\ — v(t — £1)) fixes the ST trajectory along which the particles with given 
momentum can reach the ST point x, time to is given by the initial or boundary conditions, 
v = p/p° (v = \v\), and the exponential factor accounts for the probability not to have any 
other collision from the origin point x\ till x. In the arguments of exponents n(x) is the 
particle density in the calculational frame, n(x) = J d 3 pf(x,p), o is the total scattering 
cross section. After performing integrations over d 3 x with a help of (^-functions we can 
write MBTE equation in the form: 

jfdpf = J V'il [f(t 1 ,p 3 )G(t 1 ,p 3 )f{t 2 , P4 )G{t 2 ,p A ) - f(t 1 ,p)G(t 1 ,p)f{t 2 ,p 2 )G{t 2 ,p 2 )) , (4) 
where 

V'it = \ fdh f dt 2 I D 2U WH, f(t 1 ,p)=f(t 1 ,x-v(t-t 1 ),p), (5) 

— J dt2un(t2,x—v{t—t2))v t 
ri(+ \ e 1 r<( \ f ~ ft dt 2 an{t2,x-v(t-t2))v , , 

G(ti,p) = — r , C{x,p)= dt x e Jt i . (6) 

iv[X,p) Jto 

The obvious limit in which MBTE is reduced to BTE is a completely homogeneous 
ST distribution function (i.e. no external forces, no boundaries). Another possibility is a 
thermodynamical limit, A = 1/an — > 0, when the exponential factors JHJ) will be reduced 
to ~ S(t — ti t2 ), reproducing the BTE after t x , t 2 integrations. 

For the BTE the symmetry of the invariant transition rate lead to the consequence that 
local conservation laws can be derived from the original BTE, i.e. 

d ll T ftt ' = Q and 8^ = 

were T and iV are given as momentum-integrals over the single particle PS distribution. 
Although now we have delocalized the equations, the local conservation laws can be still 
derived in the same way. Let us create a quantity ^(x) = a(x) + b(x)p%, which is 
conserved in the binary collisions 12 — > 34, i.e. + \&2 = ^3 + ^4- Now let us study the 
quantity F: 

—^yiPid»f(x,Pi) = 
Pi 
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= /^*i[/(*i,J*)G(*i,P3)/^ • (7) 

Using the symmetry of the transition rate it is easy to show [3] that 

F = J V[\ t iJ(t 1 ,p 3 )G(t 1 ,p 3 ) f(t 2 ,p 4 )G(t 2 ,p 4 ) + tf 2 (a;) - tf 3 (x) - * 4 (x)) = . (8) 

Now if we choose ^ = q, where q is a conserved charge, we obtain the charge conservation, 
and if we choose = p v k we obtain the energy-momentum conservation. 

The very essential property of the BTE is the Boltzmann H-theorem. In order to study 
the entropy 4-current, 

S»= [ ^ I fJ(x,p 1 )(log(f(x,p 1 ))-l) , 
j p 1 

we have to choose ^k{x) = log (f(x,Pk)), which is not a conserved quantity. Nevertheless, 
repeating the same steps as for eq. (|SJ) we obtain: 

^ = / V^ZJit^x-v^t - ti),j? 3 )Gf(ti,p 3 ) f(t 2 ,x-v 4 (t - t 2 ),p 4 )G(t 2 ,p 4 )x 



x 1()g ( f(x, Pl )f(x,p 2 ) \ 

\f(x,P3)f(x j p 4 )J 

The usual way of proving that Sft > does not work for MBTE because of the delocal- 
ized integral kernel, f(t x ,x-v 3 {t - t 1 ),p 3 )G(t 1 ,p 3 ) f{t 2 ,x-v 4 {t - t 2 ),p 4 )G(t 2 ,p 4 ). The 
behaviour of the entropy current in MBTE is a subject of the future studies. Neverthe- 
less, the condition of the adiabatic expansion (S 1 ^ = 0) is the same as for BTE, namely 
f(x,Pi)f(x,p 2 ) = f(x,p 3 )f(x,p 4 ). 

In conclusion, we have shown that the basic assumptions of BTE are not satisfied 
during the FO process and other very fast processes, and, thus, the description should 
be modified if used for such a modeling. We suggest the Modified Boltzmann Transport 
Equation and study some of its properties. Although the use of MBTE makes kinetical 
FO description much more complicated, some simplifications can be done [3] and the 
simplified model can be used for the qualitative understanding of the basic FO features. 
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